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Abstract. - We study the scaling of the localization length of two interacting particles in a one- 
dimensional random lattice with the single particle localization length. We obtain several regimes, 
among them one interesting weak Fock space disorder regime. In this regime we derive a weak 
logarithmic scaling law. Numerical data support the absence of any strong enhancement of the 
two particle localization length. 



Introduction. — Quantum single particle dynamics 
in one-dimensional disordered lattices with uncorrclated 
random onsite energies exhibits Anderson localization [I] . 
The asymptotic spatial decay of an eigenvector is expo- 
nential and given by Af ' ~ e~ 1 ^^ , where £X is the local- 
ization length of an eigenmode v with the eigenvalue A„, 
and the integer I counts the lattice site (see also e.g. [2]). 
The localization length is bounded from above. 

The interplay of disorder and interaction of two interact- 
ing particles (TIP) in a random one-dimensional potential 
stirred was considered by Shepelyansky (Sh94) [3J. The 
conclusion was that two particles might propagate coher- 
ently over distances much larger than the single particle 
localization length £i , if both particles are launched within 
a distance of £x from each other. Sh94 used an analogy 
between the two-particle eigenvalue problem and that of 
banded random matrices, and made an assumption about 
the scaling properties of overlap integrals which connect 
different noninteracting Fock eigenstates in the presence 
of interaction. He finally concluded that in the weak dis- 
order limit £i — > oo the two-particle localization length 
£2 will scale with £1 as £2 ^ £fU 2 , where U is the inter- 
action strength [3j. This result was further supported by 
Imry (Im95) in [4], where a Thouless-type scaling argu- 
ment was replacing the banded random matrix analogy. 
Therefore, two interacting particles were predicted to ex- 
plore a much larger space than noninteracting particles. 
Numerical calculations by Frahm et al (FR) [5] concluded 



that the scaling is probably weaker, namely £2 oc £i' 65 , 
and raised doubts about the previously assumed scaling 
properties of overlap integrals. Using a Green function 
method adapted to the problem, a new scaling relation at 
the center of the band, £2 = £1 +c£ 2 \U\/V, was derived in 
[B] (c 0.11 for bosons). In particular, this implies that 
the enhancement effect will set in for weaker interactions 
than previously predicted. Later on, it was argued that 
the enhancement effect is probably due to finite-size effects 
and it should completely vanish for an infinite system [7J. 
Simulating the time dependent Schrddingcr equation for 
two interacting particles [5j , it was argued that the dynam- 
ics is characterized by two time scales, t\ and t%, set by, 
respectively, two localization lenghts, £1 and £ 2 - Recently, 
two of us studied statistical properties of the overlap in- 
tegrals perturbatively and numerically for weak disorder 
[S] . These results contradict previous assumptions of Sh94 
and Im95 [3J[3], and if used within the previously applied 
theoretical schemes, predict a much weaker interaction in- 
duced increase of the localization length than previously 
discussed. Despite a number of studies, the problem of two 
interacting particles in a random potential remains there- 
fore a completely open problem. At the same time this 
seemingly academic case can be both addressed by cur- 
rent techniques with ultracold interacting atoms |10| , and 
is of fundamental importance for tackling the much more 
complicated case of many interacting particles in random 
potentials. 
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In the present work we first show that a nonpertur- 
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bative strong localization length enhancement can be ex- 
pected only in a regime of very weak disorder, with up- 
per bounds on the disorder strength. This regime was 
not fully accessed in previous numerical scaling studies. 
We then obtain upper bounds on the strength of the ex- 
pected enhancement effect using correct scaling proper- 
ties of overlap integrals. We then perform direct numer- 
ical measurements of the two-particle localization length 
in the perturbative and nonperturbativc regimes, by solv- 
ing the corresponding eigenvalue problem with subsequent 
averaging over many disorder realizations. Finally we for- 
mulate a set of open issues which have to be addressed in 
the future. 

Model. — We consider the Bose-Hubbard Hamilto- 
nian with disorder 



for the noninteracting problem, \q) = Y^a<v IMj v )i 
where the coefficients cp0 satisfy the eigenvalue problem 
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H — "Ho + Hint , Hint — y ' 

I 

Ho = ^ bltifvi + V + afai +1 )] , 



(1) 



which consists of non-interacting and interacting parts, Ho 
and Hi n t ■ Here a ; + and hi are standard boson creation and 
annihilation operators on a lattice site I and U measures 
the interaction strength. The random on-site energies e; 
are chosen uniformly from the interval [-W/2, W/2], with 
W and V denoting the disorder and hopping strengths, 
respectively. 

One particle. In this case the interaction term does 
not contribute. We use the basis \l) = af\0) with I = 
1, . . . , N (N is the number of lattice sites) . The eigen- 
states (also called single particle normal modes (NM)) 



Y^f A i u '\l) are defined through the eigenvectors 



A\ ~ e with the eigenvalue problem 

A ^)=e^ ) +nA ( : ) i+A ( - ) i)- 



(2) 



The eigenvalues -2V - W/2 < A < 2V + W/2 fill a band 
with a width A x = W + W. The most extended NMs 
correspond to the band center A = with localization 
length 

£i(A = 0, W) « 100{V 2 /W 2 ), (3) 

in the limit of weak disorder W/V < 4 [2J. The average 
volume L which an eigenstate occupies has been estimated 
to be about L ps 3£ x for weak disorder [9]. 

Two particles. For U = we construct orthonormal- 
ized two particle eigenstates as product states of single 
particle eigenstates in a corresponding Fock space 



\fj,,v> fj) 



V 1 + <w 



H \fi,v) = (A M + A„)|/i,i/). (4) 



Then, we expand the eigenstates |g) of the interacting 
particle problem, H\q) = A g |g), in systems of eigenstates 



(5) 



u' ,v' 



Here A M „ = X u + A„ and therefore the noninteracting case 
U = yields a band with width A 2 = 2A X . The coeffi- 
cients Ifi" are connected with the overlap integrals 



Tl*'"' — 4 M A v A 11 ' A v ' 



as follows: 



If 1 



(6) 
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The interacting case yields a single band for U < A2, but 
two bands separated by a gap for U > A 2 . Indeed, in the 
latter case two-particle bound states are renormalized out 
of the main band, and are mainly consisting of two par- 
ticles occupying the same site Therefore, remaining 
band is due to states where the two particles can be any- 
where but not on the same site. This is simply the limit 
of two noninteracting spinlcss fcrmions. The localization 
length of these two noninteracting fcrmions is of the same 
order as the single particle localization length. The lo- 
calization length in the bound state band is even smaller, 
since the effective disorder strength in this band becomes 
2W, but the effective hopping is strongly suppressed. 

For numerical purposes we expand the two particle 
eigenstates \q) in the local basis 
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where C^L = (l,m\q) arc the normalized eigenvectors. 
They satisfy 



_ \ - A m -flj + /jj A m r q ) 

w — m~i — / -I 1 r — ' i,m 
L rr v 1 + dim ■ \/l + o av 

rn,l<rn v V 



(9) 



We will numerically compute the probability density 
function (PDF) of the number of particles in direct space 
pi = (q\afai\q) /2, which is given by 



Pi 



(./) 



/ N N \ 

\ E 4f + E 4S 2 

\k,l<k m,l>7n J 



(10) 



On scales. — Since a single particle eigenstate occu- 
pies a volume L, there are of the order of L 2 two parti- 
cle eigenstates which are residing in the same volume for 
U = 0. The overlap integrals built among these L 2 Fock 
states are nonzero (more precisely not exponentially weak) 
and define the connectivity in the Fock space for nonzero 
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U . The average eigenenergy spacing d of these connected 
Fock states is d = A2/L 2 . It therefore defines an effec- 
tive energy mismatch - i.e. an effective disorder strength 
W = d - in the Fock space. The effective hopping strength 
follows from ([S]) and is given by V — 2U (I). Here (/) is an 
average overlap integral among all connected Fock states 
0. 

In analogy with eq. (|3|) we can therefore obtain a local- 
ization length in Fock space for weak Fock space disorder 
W < 4V, which in real space is a measure in units of the 
single particle localization length: 

!«100|^ 400^</) 2 i 4 . (11) 

For strong Fock space disorder W ^> V the volume 

and two interacting particles are localized in the same way, 

therefore £2 ~ £1 in this case. 

Bounds on the weak Fock space disorder regime. Let 
us now address the question whether we can enter the weak 
Fock space disorder regime for strong single particle disor- 
der W>K This seems possible at a first glance since we 
can increase the value of V by increasing U. However, in 
this limit (I) ~ V 2 /W 2 . Therefore the needed interaction 
strength is U ~ W 3 /V 2 , since A2 ~ W. But an increase of 
the interaction strength beyond the band width A2 leads 
to the separation of the energy spectrum into two bands - 
a bound state band with strongly localized particle pairs 
|llj , and a noninteracting spinless fcrmion band which has 
no localization length increase as compared to the single 
particle case. The two conditions U < W and U > W 3 /V 2 
imply that W < V is needed, which means that the single 
particle case must be in the regime of weak localization. 
Therefore U < V is an upper bound for entering the weak 
Fock space disorder regime. 

Lowering U further we will however again leave this 
regime and enter the pcrturbative one, which is again char- 
acterized by strong disorder in Fock space. Indeed, the 
energy renormalization of a given Fock state follows from 
lp)]) and is given by 2UIq where 1$ is an average overlap 
integral of a Fock state with itself. Due to orthonormality 
of the single particle eigenfunctions it follows Iq w 1/L. 
The perturbative regime holds as long as UIo < d. Inside 
the perturbative regime a Fock state is still a good ap- 
proximation to an exact eigenstate, and therefore the two 
particle localization length is of the order of the single par- 
ticle one. Therefore, the nonperturbative weak Fock space 
disorder regime is accessed for A2/L < U < V. 

For any practical purposes we seek a strong enough in- 
teraction strength U, and this requires U ~V and W < V. 
In order to obtain any relevant scaling results upon vari- 
ation of W one needs therefore to lower W significantly 
further such that W < V. 

Overlap integrals revisited. Sh94 and Im95 estimated 
the average overlap integral (I) si ~ L~ 3 / 2 pHH] inside 
the weak Fock space disorder regime. This result is ob- 
tained in the following way. A single particle eigenstate 



occupies a volume L> 1. Due to normalization it follows 
\A\ u) \ ~ L- 1 ' 2 . The crucial point was to assume that 
all terms inside one localization volume in the sum ([5]) 
have uncorrclatcd signs. This leads to the above estimate. 
However, in the limit of weak disorder and large localiza- 
tion length, the single particle eigenvectors inside a local- 
ization volume will appear similar to plane waves, with 
appreciable phase correlations between different sites, and 
also between different eigenstates. Some numerical studies 
by Romer et al (R99) Q~2] even concluded that (I)r ~ L~ 2 . 
This result essentially corresponds to the assumption that 
the eigenvectors are exact plane wave states inside a local- 
ization volume. It is this small difference in the exponent 
which separates a possible existing strong enhancement of 
the localization length from no effect at all. 

In a recent work two of us performed a perturbation 
approach at the weak disorder limit and obtained that 
strong phase correlations will certainly modify the predic- 
tion of Sh94,Im95. At the same time corrections to the 
result of R99 arc significant. As a final result we obtain 
CO SJ ~ — ln(L)L -2 [Hj - logarithmic corrections to the 
prediction of Romer et al. It is well-known that logarith- 
mic corrections are rather resistent to numerical verifica- 
tions, if no special trick or technique is used. Therefore, 
our numerical tests in a limited interval of W lead only to 
the clear result that the prediction of Sh94,Im95 is incor- 
rect, and if (I) ~ 1/L Z is assumed, then z ps 1.7. They 
were not sensitive to distinguish between this power law 
and a possible asymptotic (I) ~ — In (L)L~ 2 logarithmic 
law. 

Scaling of the localization length. Combining the 
above predictions on the overlap integral scaling and the 
localization length scaling (fTTj) we arrive at the following 
results in the weak Fock space disorder regime. Here we 
set A 2 = 8V, take W < W such that © holds. Then 
Sh94 and Im95 predict ~ {U/V) 2 ^i as derived using 

different methods in the original papers 013]. According 
to R99 the whole effect is simply £2/^1 ~ (U /V) 2 , i.e. no 
enhancement at all. Finally, our analytical estimate for 
the overlap integrals yields 

Note that the numerically estimated overlap integral de- 
pendence on L results in £2/^1 ~ {U/V) 2 ^- 6 . 

Numerical technique. — We estimate the largest av- 
erage localization length £2 of the probability density func- 
tion pi ~ e ~ 21 ^ 2 [see Eq. (fTO)) ] using the following proce- 
dure (the prefactor 2 in the exponent takes care of the 
fact that densitites instead of wave functions are fitted). 
For a given realization we solve the eigenvalue problem 
and choose only those modes df} yl which satisfy to the 
following selection rules: 
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Fig. 1: (a) Smoothed average probability distribution function 
(pi) versus lattice site I in lin-log scale for W = 2, U = [(g), 
green curve] and W = 2, U = 0.2 [(o), orange curve]; (b) the 
corresponding quantity a (see text) versus I, with a zoom of the 
interval with saturated values of a (inset) ; (c) the two-particle 
localization length £2 versus W for the noninteracting case, 
U=0 (red circles). Blue solid line: £1 = 100/ W 2 . Dashed lines: 
maximal admissible error of 10% from the analytical formula. 
Gray area corresponds to the admissible values. 

• the center of masses 

N N 

k - E lc £ . ™*= E mC{ £ ( 13 ) 

m,/<m m,l<~m 

satisfy to the inequalities \l q — N/2\ < £1 \fh q —N/2\ < 
£1 (( is of the order of the corresponding average lo- 
calization length for a single particle problem). Thus, 
we take into account only those modes for which the 
two particles reside in the same localization volume; 

• the eigenvalues are near the bandwidth center. We 




Fig. 2: The ratio &/£i versus interaction constant U for dif- 
ferent values of disorder W = 2, 2.5, 3, 4 (from top to bottom) 

assume that similar to the case of a single particle 
problem the most extended modes are with X q 0; 

• we project £; onto the modes of the one-particle 
problem, calculate the amplitudes <^ M „ in accordance 
with Eq. ^ and find the mode \/j,o,vq) with the 
largest amplitude, max^^,^. Such a method al- 
lows us to identify the Fock state |/xo, vq) which domi- 
nates all others. We then request that the eigenvalues 
A M0 and A„ are close to the bandwidth center. Thus, 
we exclude possible cases when X q is close to the band 
center, but and are located at the two oppo- 
site band edges. 

Having selected the modes C/j q ^ , we compute their proba- 
bility density functions pi according to Eq. (fTU)) and shift 
them such that their new center of mass are located at 
the center of a chain, N/2. Then, we compute loga- 
rithms of the PDFs, in(p/) and perform a statistical aver- 
age of the PDFs over many disorder realizations as (pi) = 
cxp[(ln(pj))]. Finally, using a local regression smoothing 
technique, we obtain smooth functional dependencies of 
(pi) and calculate the quantity a = 2 ■ |d(ln(p;))/d/| _1 . In 
the limit of large I, a (I) should saturate at the average two 
particle localization length £2- 

Numerical results. — The dimension of the Hilbert 
space p grows rapidly (~ TV 2 ) with the size of a chain, so 
that the maximal reachable size used in numerical compu- 
tations, N max = 234. Thus, we inevitably face finite size 
effects for weak disorder. We start with the noninteracting 
case U — for which £2 must be exactly equal to £1. We 
estimate the minimal value for the strength of disorder, 
respectively, maximal localization length, £2, at which an 
error (caused by finite-size effects) is less than 10% (which 
is the maximal error we admit). We assume that this error 
depends only on the magnitude of £2 but not on the inter- 
action strength U. Thus, the largest tolerable values for 
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Fig. 3: The two-particle localization length £2 versus one par- 
ticle localization length £1 for U = 1.5, 1.0, 0.5, 0.2 (from top to 
bottom) in log-log scale. Dashed straight lines are power laws 
£f with the exponents a = 1.4 (upper line) and a = 1.3 (lower 
line). The size of a chain is N = 234. 
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Fig. 4: log(|ci jm | 2 ) versus I and m averaged over time and 
5000 disorder realizations. For better visualization results are 
unfolded from the irreducible triangle shaped state space onto 
a square with c; jm = c m j for m < I. The strength of disorder 
W = 2.5 and the interaction constants U = (a) and U = 2 
(b) . Particles are initially located on the same site at the center 
of a chain with N = 170 sites. 



£2 found for the nonintcracting case are also assumed to 
be the limiting values for the interacting case. For U = 
the lower curve in Fig.flja) presents a smooth dependence 
of (pi) on I. The corresponding quantity a (lower curve in 
Fig- Hlb)) saturates at large distances. The obtained local- 
ization length £2 is finally shown in Fig. [TJc) and agrees 
well with the theoretical prediction, however systematic 
deviations accumulate for weak disorder. A recalculation 
of the same quantities for U = 0.2 in Fig. [IJa,b) shows 
that the method appears to be applicable to the inter- 
acting case as well. Finite size effects blurr our results 
substantially if £2 > 40. 

Let us discuss our results for nonzero interaction. The 
ratio £2/^1 grows with increasing interaction constant U , 
as shown for different values of W in Fig. [2] This growth 
is stronger, the weaker the disorder strength is. For our 
data, the ratio did not substantially exceed the value 2. 
However, it seems plausible that for FT < 2 (which is not 
treatable with our current technique), stronger enhance- 
ment effects could be observed. 

The central result is plotted in Fig. [3] Here we plot 
£2 versus £1 on log-log scales. We try to fit data for a 
fixed value of U and different values of W using power law 
estimates. Both £1 and £2 vary less than an order of mag- 
nitude, while a safe power law fit needs at least two orders 
of magnitude variations on each variable. Nevertheless we 
bound the obtained variations with two lines £2 ~ £l' 3 
and £2 ~ £i' 4 - Such a scaling is much weaker than the any 
of the above predicted power laws. It is possile that we 
observe the onset of the logarithmic scaling obtained from 
perturbation theory (p~2j) . 

Averaged evolution of two particles. In order to vi- 
sualize the effect of interaction on the localization of two 



particles, we solve the time dependent Schrodingcr equa- 
tion idt\^(t)) = H\^(t)). We expand |#(i)) in terms of 
the orthonormal states \l,m) (I < m) as 

N 

!*(*))= E %m{t)\l,m)- (14) 

rn , I < rn 

where the coefficients C\ q ^ n 

Q, ro (o=E^4t e " iV - ( 15 ) 

Here ip q are the amplitudes of NMs related with the initial 
amplitudes q, m (0) = (I, m\^(0)) of the two-particle states 
as 

N 

E Cl, m (Q)4%. (16) 

m,/<m 

We launch two particles on the same site, Iq = mo, or 
adjacent sites, Iq = mo — 1, such that the initial amplitude 
c;, m (0) = 5Li a o~ m ,m - We calculate then the averaged in 
time square amplitude (|c; iTO | 2 ) t [see Eqs. (|15|16)) ]. which 
is given by 

(\c hm \\ Mm £^Ldt = £ |^| 2 4t 2 . (17) 

9=1 

We further average (|Q, m | 2 )t over 5000 disorder realiza- 
tions. In addition we perform an averaging with re- 
spect to initial conditions, by keeping the same disor- 
der potential, and taking different neighboring sites as 
an initial location of the particles. Finally, we com- 
pute the average probability density function (pi) using 

Pi = 2 (Sw<fc l Qfc l 2 +El,i> ra \c m i\ 2 y 
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Note that the averaged in time two-particle wavefunc- 
tion \ci m \ 2 for a single disorder realization has many spots 
at different locations due to resonances. This feature is 
smeared out, once the averaging with respect to disorder 
realizations is performed as is seen in Fig. BJa),(b). For 
the noninteracting case the obtained distribution is elon- 
gated along the main axes. This happens because the two 
particles are not correlated, and it is much more proba- 
ble for them to occupy different space regions. However 
for U = 2 the distribution is elongated along the diago- 
nal. This implies that the two particles are exploring more 
states when being close to each other. 

Summary. — In summary, wc discussed the possible 
regimes of two interacting particles in a random poten- 
tial. The most interesting case of a weak Fock space disor- 
der regime was analyzed, and scaling laws were discussed. 
These results, as well as the numerical data presented as 
well, show that the localization length enhancement effect 
is much weaker than previously assumed. Further numer- 
ical studies are needed in order to substantiate these re- 
sults. However the current techniques are not of use for 
weaker disorder strength. Therefore new computational 
approaches arc needed in order to reach disorder values as 
low as W = 0.1, which may be enough to test the pre- 
dicted weak logarithmic scaling. 

Acknowledgement. — The authors wish to thank I. 
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